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Abstract 

The confluent second-order supersymmetric quantum mechanics, with factorization 
energies e±, e<i tending to a single e- value, is studied. We show that the Wronskian 
formula remains valid if generalized eigenfunctions are taken as seed solutions. The 
confluent algorithm is used to generate SUSY partners of the Coulomb potential. 
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1 Introduction. 

The supersymmetric quantum mechanics (SUSY QM) represents a powerful tool for the 
spectral design in quantum theory. It has proved successful to construct solvable potentials 
for which the spectral information can be anallytically determined [1-5]. The simplest version 
is the first-order SUSY QM which, however, does not allow to modify levels different from 
the ground state without creating singularities in the SUSY partner potential [3]. As an 
alternative to surpass this problem, the second-order SUSY QM can be used [6-13]. Indeed, 
through this technique one can embed two levels e%, e 2 between two neighbouring energies 
of the initial Hamiltonian [5,7]. It is possible as well to create single levels above E [13] 
and to generate complex potentials with either purely real spectra or having some complex 
'energies' [12]. 

For ei 7^ £2 the modification to the potential induced by the second-order SUSY QM 
involves the Wronskian of the associated seed Schrodinger solutions u\, U2- In the confluent 
case ei = e 2 = e, however, the key function w(x) inducing the change depends just of one 
seed solution u and it is not clear that the treatment based on the Wronskian remains valid. 
In this paper we will prove that result by identifying appropriate (generalized) eigenfunctions 



for which the Wronskian coincides with w(x). This implies that a unified technique valid in 
the non-confluent and in the confluent case is available, we just have to identify the right 
seed solutions. 

In the next section we will address the second-order SUSY QM and its classification 
scheme. We will study in some detail the confluent case, the restrictions onto the seed 
solution to ensure the regularity of the potentials difference and the discussion about the 
Wronskian. Then, we will generate confluent second-order SUSY partners of the Coulomb 
potential (some of which are new). We will end the paper with our conclusions. 



2 Second-order supersymmetric quantum mechanics. 

The second-order SUSY QM consists in the following realization of the standard SUSY 
algebra with two generators: 

{Qj, Qk} = S jk H ss , [H„, Qj] = 0, j, k = 1, 2 (1) 



A = 17i + i {x) Tx +l(x) < Sa = ah < 5 > 

The intertwining relationship (j5J) leads to the following set of equations: 

~2 ^ + v[^-V-V + dj+c = (6) 

2 / 

7 = d-y + l.-l, V = V + 2 V ' (7) 

Given V(x), c and d, the new potential V(x) is determined by the solutions r\ to the non- 
linear equation (JHJ). To find them, let us use the Ansatz: 

rj[{x) = r] 2 (x) + 2(3{x)r]{x) - 2£(x) (8) 

which leads to £ 2 = c and the following Riccati equation: 

p + (3 2 = V - e, e = d + f (9) 

We can work instead with the equivalent Schrodinger equation (fl = u'/u) 

-u" + Vu = eu (10) 

Thus, the solutions r] can be classified according to the sign of c. 



2 



3 Classification of second-order SUSY transformations. 



3.1 The non-confluent cases with c 7^ 0. 

Here e\ = d + t/c, e 2 = d — y/c, €\ 7^ e 2 ; this includes the real case with c > 0, ei, €2 G K and 
the complex one with c < 0, t\ G C, e 2 — e*. From (jHJ) two equations for 77 are obtained: 

V ' = V 2 + 2(3 lV - (e 1 - e 2 ), V = r/ 2 + 2(3 2 V + (ei - e 2 ) 

which leads to: 

where g) = /(/ — gf is the Wronskian of / and g. To avoid singularities in rj(x) 

the Wronskian has to be nodeless. The spectrum of H, Sp(if), will depend as well on the 
normalizability of the two eigenfunctions ij) e . of H with eigenvalue 6j in the Kernel of A', 

with explicit expressions given by ipa °^ U2/W(ui,u 2 ), ip e2 °^ Ui/W(ui, tt 2 ). 

Some spectral design possibilities are worth to be mentioned (for a detailed treatment see 
[5]). i) Two levels e±, e 2 can be created below the ground state energy of H, namely, e\ < e 2 < 
Eo, Sp(i^) = {ei, e 2 , E n , n = 0, 1, 2, • • • } [8]. A pair of levels can be placed between two 
neighbouring energies of H, i.e., < e x < e 2 < E i+1 , Sp(if) = {^ , • • • , Ei, e x , e 2 , • ■ • } 
[7]. Hi) Two neighbouring energies of H can be deleted, namely, e\ = E iy e 2 = E i+1 , 
Sp(H) = {E , • • • , Ei_i, E i+2 , ■ ■ ■} [5]. iuj Some complex energies can be manufactured (the 
new Hamiltonians are non-hermitian) [12]. 

Cases ii), in), iv) run against the dominant idea that in SUSY QM the new levels are 
always real and below the ground state energy of the initial Hamiltonian, improving thus our 
spectral design possibilities [4,5]. 



3.2 The confluent case with c = 0. 

We have now that e = e\ — e 2 , and the Ansatz (JHJ) just provides [10, 13]: 

7/ = rj 1 + 2(5r] 

This Bernoulli equation has a general solution given by 

g 2 / p(x)dx j 

V{X) = Wo -fe*f^d X = ~T X lnHx)] (12) 
where, up to an unimportant constant factor, 

w(x) =w - u 2 (y)dy (13) 



Since equations ()lipi2j) look similar, perhaps w(x) has to do with a Wronskian. To see that, 
suppose there is a function v related with the given u through: 

(H-e)v = u, {H-e)u = Q (14) 
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i.e., u and v are rank 1 and rank 2 generalized eigenfunctions of H. This fact was noticed 
in [4], but a further study exploring the link with the Wronskian formula is needed. Let us 
find now v from the differential equation t)14j) : 

e=u ( t+ /3£H (i5) 

Up to a constant factor, it turns out that w(x) = W(u, v). 

Similarly as for t\ ^ 62, the regular confluent transformation will be produced by a 
nodeless w(x). If the domain of x is M, it is sufficient to use solutions u(x) vanishing either 
when x — ► —00 or when x —>■ 00 [13]. However, if the x-domain is the positive semi-axis, as 
in the Coulomb problem which we will address later, the simplest choice is to take: 

limw(x) = or lim u(x) = (16) 

a;— »0 x^oo 

For e G" Sp(-ff) one of the requirements (|T5j) can be achieved, but not both simultaneously. For 
e G Sp(-ff) the two conditons (jTfij) are automatically satisfied. In both normalized 
eigenfunction of H with eigenvalue e in the Kernel of can be found, ip e oc u/w. In 
particular, for e> E equation (fTB^) can be fulfilled, i.e., one level can be created above So- 
Let us remind that the excited state solutions obey (jTUjl and thus they are appropriate to 
generate new potentials through the confluent SUSY algorithm. 



4 Confluent SUSY partners of the Coulomb potential. 

Let us apply the confluent algorithm to the Coulomb problem which, after separating the 
angular variables and taking H — e = m — 1, leads to a Schrodinger Hamiltonian as given in 
(J3|, with an effective potential 

„ (r) _-? + 3£+ii (17, 

where I = 0, 1, ... . The spectrum of the corresponding Hamiltonian reads Sp(H) = {E n = 
— 1/n 2 , n — £ + 1,1 + 2, ... }. First we find the general solution to (fTUj) with the V(r) of (fTTj) 
and an arbitrary factorization energy e < 0; second, the solution with the right behaviour 
is chosen (see Eq. (116(0 . After implementing the procedure, the solution vanishing at the 
origin reads: 



u\r) 



(-f^+Tb) (2r^Y +1 e W=i 1 Fx U + 1 - 21 + 2; 2r^) (11 
' 1 W(2< + 2)p V v ; y ^> ! v y v 



T(z) and iFi(a,b;z) being the Gamma and the confluent hypergeometric functions respec- 
tively. Notice that u(r) becomes the normalized eigenfunction of H for e = E n . A straight- 
forward calculation leads now to our key w-function: 



■B{ l + l + -^.C + l +m -^-)(2r x f^ 

w(r) = w — 2 



111: 



1o 2(2£ + m + 3)(2£ + 1)1 ml B -1,1+ 1 - -j= 



x 2 F 2 ( 2i + m + 3, £ + 1 + — =; 2£ + m + 4, 2£ + 2; -2rv^ J (19) 
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where B(x,y), 2-^2(^1, 02! &i> &2! z) are the Beta and a generalized hypergeometric function. 
For e = E n = — 1/n 2 this infinite series truncates: 

w(r) = wq — 



m=0 



(2£ + m + 3)(2£ + m + 1) 



2 F 2 (2l + m + 3, n + £ + 1; 21 + m + 4, 21 + 2; -f ) 
X 2n(2£ + l)!m! J B(n-£-m,2^ + m + l) ^ ^ 

The wo-domain where w{r) is nodeless reads (— 00, 0] for e ^ E n and (— 00, 0] U [l,oo) for 
e = E n . 

The SUSY partner potentials take the form: 

k W = J + 5M + ^ + ^ (21) 

r H w(r) w 2 (r) 

In (|21j) w(r), w(r) are given by (|18j) and (jl°|20j) while u'(r) arises from (|T8|) . Some examples, 
induced by the w(r) of (ffilj) . deserve an explicit discussion. 

4.1 The case with £ = and n = 1. 

Here the SUSY partner potentials of V^(r) = — 2/r are given by: 

y (r) _ 2 16 r [-l-r+( Wo -l)(r-l)e 2 -] 
r [1 + 2r + 2r 2 + (w - l)e 2 *f 

The potentials V(r), V(r) are isospectral for w; G (— 00, 0) U (1, 00) and they differ in the 
ground state for u> = 0, 1 because E\ = — 1 is missing of H in the last case. The potentials 
()22|) coincide with a family derived long ago through the factorization method, we just take 
w = 47x in (J22D to get equation (3.1) of [14] (see also [12,15,16]). 

4.2 The case with £ = and n = 2. 

The SUSY partner potentials of V(r) = —2/r become now 

~ 2 8r(r-2)[-8 + 4r + 6r 2 + 2r 3 + r 4 -2(w -l)(4-6r + r 2 )e r ] 

r [8 + 8r + 4r 2 + r 4 + 8(w - l)e r ] 2 

Once again, V(r) and V(r) are isospectral for w G (—00, 0) U (1, 00), and for w = 0,1 their 
spectra differ because the level E 2 = —1/4 is missing of H. As far as we know, the potentials 
(J23|) have not been derived previously. 

Other families of exactly solvable potentials can be written explicitly for different values 
of n and £. We illustrate one of them in figure 1 for n = 4, £ = 1 and wq = —0.1. For 
comparison, the initial potential is also shown. 
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Figure 1: Isospectral SUSY partner potentials V(x) (black curve) and V(r) = — 2/r (gray 
curve) generated through (fTHj). ()20I21|) with n — 4, £ = 1, Wq = —0.1. 

5 Conclusions 

Contrasting with the first-order SUSY QM, in which we modify just the ground state energy, 
the confluent algorithm allows the embedding of levels at any positions on the energy axis. 
We have shown that the Wronskian formula is still valid in the confluent case. We conclude 
that the second-order SUSY QM represents a poweful tool in which the right choice of the 
seed solutions allows us to generate the new potential in a simple way. 
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